Indonesian Journal of Electrical Engineering and Computer Science 
Vol. 27, No. 3, September 2022, pp. 1151~1161 
ISSN: 2502-4752, DOI: 10.1159 1/ijeecs.v27.i3.pp1151-1161 Oo 1151 


Improved security and stability of grid connected the wind 
energy conversion system by unified power flow controller 


I Made Wartana!, Ni Putu Agustini!, Sasidharan Sreedharan? 
'Department of Electrical Engineering, National Institute of Technology Malang, Malang, Indonesia 
*Department of Electrical and Electronics Engineering, MES College of Engineering, Kuttippuram, India 


Article Info 


ABSTRACT 


Article history: 


Received Nov 17, 2021 
Revised May 25, 2022 
Accepted Jun 16, 2022 


Keywords: 


Doubly-fed induction generator 
Max. LBS 

Min. Pioss 

Small-signal stability 

Unified power flow controller 


The stability and security improvements of the grid-connected to the wind 
energy conversion system (WECS) can be made by optimizing the 
placement of a flexible alternating current transmission system (FACTS). 
This study discusses the optimal placement of one type of WECS, namely 
the doubly-fed induction generator (DFIG) with a series and a shunt-FACTS 
control device called unified power flow controller (UPFC). The DFIG and 
UPFC connected grid dynamic performance improvement with a maximum 
load bus system scenario. The optimal placement of DFIG and UPFC on the 
grid is formulated as a multi-objective problem, namely maximizing load 
bus system (Max. LBS) while minimizing active power loss (Min. Pioss) by 
pleasing numerous security and stability constraints. The non-dominated 
sorting genetic algorithm II (NSGA-II) have been utilized to settle this 
opposed bi-objective enhancement problem. The validity of the suggested 
method was examined on a modified IEEE 14-bus and a utilitarian examine 


system connected to DFIG with UPFC in power system analysis toolbox 
(PSAT) software. The optimal placement of DFIG and UPFC on the grid has 
increased the system's dynamic performance, with all the specified particular 
constraints being encountered. 
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1. INTRODUCTION 

Electrical energy plants distraibute electrical power to meet consumer needs by prioritizing safety, 
reliability, and stability [1], [2]. Recently, conventional energy generation technology, which mainly uses 
coal, petroleum, and nuclear fuel, has experienced a setback due to the depleting availability of fuel oil, 
economic awareness, concern for environmental issues, and safety measures. Therefore, the development of 
efficient renewable energy conversion and utilization techniques is the main focus in the future. One of the 
efficient conversion techniques is wind energy because it utilizes environmentally friendly and pollution-free 
energy sources [3]. 

Therefore, the number of wind energy conversion systems (WECS) connected to the system has 
increased rapidly. Some problems arise due to the integration of the WECS to the grid, such as the double-fed 
induction generator (DFIG), and need detailed analysis in solving these problems [4]. Integrated DFIG into 
the grid requires installing a shunt controller to support voltage requirements and a series controller to control 
line flow, increasing the system's load capacity. The integrated power flow controller (UPFC) can 
simultaneously control line flow and voltage [5], [6]. Therefore, to solve the problem of maximizing the load 
bus system (Max. LBS) while reducing the active power loss (Min. Pioss), it is necessary to consider various 
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security and system stability constraints [7]. The optimal placement and settings of UPFC to minimize active 
power loss in an optimally integrated DFIG in the grid to maximize the load bus system need to be 
determined. 

Numerous challenges arise due to the large-scale integration of renewable energy resources into the 
network to improve system performance [8], [9]. Sisay and Jately [10], to enhance the dynamic performance 
of grid-connected wind farms and improve the system's transmission capability, a UPFC is interconnected at 
the voltage at the point of mutual pairing. The use of modern UPFC integrated into the network to improve 
the dynamic performance of the wind farm was proposed in [11] and [12] to reduce power quality problems 
by increasing the transmission capability of the dynamic response system of the ac network. Liu et al. [13], 
the stability of the integrated grid system with DFIG and direct drive permanent magnet synchronous 
generator and UPFC was analyzed by comparing the structure, principles, and mathematical simulation 
models of the two wind turbines. The ant-colony optimization (ACO) and adaptive neuro fuzzy inference 
system (ANFIS) algorithms are proposed in [14] for the WECS and UPFC optimal reactive power delivery 
(ORPD) techniques in grid-connected power systems. Integrated power system using DFIG with three-phase 
fault using UPFC [15], whereas, Vig and Surjan [14] proposes a high-precision new grid-connected three- 
phase inverter system with a sub-inverter using a combination type UPFC structure consisting of a central 
inverter, a sub-inverter, and a transformer. The stability of DFIG and UPFC integrated into the power system 
on the critical clearing fault time in a multi-machine power system has been investigated in work [16]. 

This paper deliberates integrating WECS type DFIG, and flexible alternating current transmission 
system (FACTS) devices type UPFC into the network to improve network security and stability. The problem 
is solved by multi-objective operational control, including max. LBS strategy and min. Poss by encounter 
several line stability and security constraints such as small-signal stability, generation limit, voltage limit, and 
the loading limit. To achieve this dynamic goal, simultaneously, the most versatile series and shunt-FACTS 
controller, UPFC, have been implemented. The optimal location and sizing of the DFIG and UPFC into the 
grid is formulated as a multi-objective optimization problem implemented on a modification of the IEEE 14 
bus standard test system and the practical Indonesian Java-Bali 24-bus system. 


2. RESEARCH METHOD 
2.1. WECS modelling 

The most commonly used wind energy conversion system (WECS) nowadays is DFIG since it has 
some advantages, working at despicable wind speed synchronous speed and regulating output power and 
voltage [17]. It is composed of an induction generator and a winding rotor with two converters was to inject a 
universal voltage into the rotor circuit at a phase shift and a varying frequency. Both converters can control in 
two means, either alternating current/direct current (AC/DC) (forward) or DC/AC (reverse) in this system. 
The rotor-connected converter and the grid-connected converter are called rotor converter and line converter, 
respectively. The DFIG provides flexible operating and robust control, being power can flow in both 
directions as per system conditions. Figure 1 shows a DFIG component model that adequately captures 
electromechanical dynamics described in detail [18]. As shown in (1) formulates the generator's active and 
reactive power production depending on the stator and converter currents and the stator and converter 
voltages, respectively. 

Ph = Ysa isa + Vs.q isq + Ve alca + F (1) 


qh = Vs,qls,a = Vs,alsq + Veqtea — Vea leg 


If expression (1) is rewritten, it can be a function of the stator and rotor currents being is,a +jisq and ira +irg, 
and the stator and rotor voltages Vs,a + jVs,q and Via + jVr,q, respectively. 
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Figure 1. DFIG component modeling 
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2.2. UPFC modeling 

This type of series and shunt FACTS controller is designed for precise timing control of the power 
system and dynamic compensation for undesired conditions. It can control the power flow of the system by 
adjusting the varied factors that influence it. Because of its ability to generate a voltage at any phase angle 
and independently regulate the real power flow regardless of reactive power, the UPFC is the most 
commonly used. However, it is relatively more expensive among various FACTS device compensation 
techniques [5]. The UPFC consists of a static synchronous compensator (shunt voltage source inverter (VSI)) 
and a series static synchronous compensator (VSI Series) coupled together via a typical DC link capacitor, 
providing a two-way real power exchange between the shunt and the series section as shown in Figure 2 [19], 
[20]. By adding shunt power and serial power insertion on buses i and j, the UPFC’s basic structure and 
mathematical model can be formulated as shown in Figure 2(a) and Figure 2(b), respectively. While 
equations (2) and (3) state the equivalent power injection element [11]. How quickly and accurately the 
reference signal is generated will determine the performance of the UPFC where a suitable dc-link current 
regulator is used to establish the actual reference signal after eliminating the fluctuating signal. 


P, upre = 0,02rbseV? siny — 1. —2rbseV;V; sin( 8; — 0; + i (2) 
P; upre = TbseViVj sin( 6; — 8; +y) 
Qi upre = —rbseV? cosy } GB) 
Qj.uprc = ThseViVj cos( 6; — 6; +y) 
V, V: I bus j 
Nar 
Vi v; 
] 
Shunt Series 
= VSI VSI 
(a) (b) 


Figure 2. UPFC model in power system; (a) basic structure and (b) mathematical model 


2.3. Problem formulation 

This study elaborates conflicting bi-objective optimization on maximizing load bus system (Max. 
LBS) with minimizing active power loss (Min. Pioss) by pleasing numerous security and stability constraints 
to verify the method proposed. The optimization problem is conceivably concurrently determined by non- 
dominated sorting genetic algorithm IT (NSGA-II), which formulates as f(x,u) in (4) and (5) [21], with x and 
u represented dependent and control variables, respectively. 


Max. f (x,u) = [fi (x,u), -f2 (x, w)] (4) 
gi(x, u) = 0i = 1,...... ,m 
St ee. u) <O0j=1,...... yn (5) 


where the bi-objective functions are represented by f; and f2, with g; being the i equality, and hj is the 
j inequality constraints, respectively. 


2.4. Maximizing load bus system (Max. LBS) 

A loaded bus increase scenario À is implemented as the primary bi-objective function in this study 
by Max. LBS, as shown in (6) and (7) based on the system load parameter VL, viz the aggregate of the 
thermal (OLL;) plus bus violation limit factors (BVVj) discussed in detail in [21]. 


Max. f,(x, u, A) (6) 


st. VL = X; OLL; + XE, BVV; (7) 
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2.5. Minimizing active power loss (Min. Pioss) 

The second objective function is Min. Pioss the transmission line as shown in (8) [22], [23]. 
The voltages on bus i and bus j of the k” line are expressed as V6; and V6), respectively. Whereas N; is the 
number of transmission lines, and gi; is the the k” conductance. 


Min. fo(X,U) = Pross = Legis gen V2 + VP — 2- V; + V; + cos 0;;] (8) 


2.6. Equality and inequality constraints 
2.6.1. Equality constraints 

The equality constraints expressed in (9) are typical load flow equations where N, is the number of 
buses. The active and reactive power loads are Pz; and Qz; respectively. Moreover Gj, By and 0 are the 
admittance and angle of admittance. 


N å 
E e n) (9) 
, = L,4;...... IND 
Qe; = Qi + V; < Di Y . (Gij “SIR ij + Bij * COS 0;;j); 


2.6.2. Inequality constraints 

The inequality constraints in this study are active Pci and reactive power generators Qg;, voltage Vi, 
and phase angle 6;(10). The parameter settings of UPFC: series voltage vs and shunt current isy as shown in 
(11). Furthermore, Pj in (12) is transmission loading. 


min max 
PEs Pe Pe 


l 


orhe Qe, < gre 


l i=1,2,....,m (10) 
wen < V; < Ve; 

ae < v < vg? an 
ish’ < isu < is? 
[Py] < Py ij = toed N (12) 


2.7. Stability constraints 
2.7.1. Small-signal stability 

As shown in (13) described the state-space model of the power system intended for small signal 
stability analysis that formulated as a differential equations and algebraic equations (DAE) [24]. 


zeren m 


0 = g(x,y) 


Where x is a vector of state variables and y is an algebraic variable, respectively. By linearizing the 
state-space model using the Taylor series expansion at a stable operating point (xo, yo), as shown in (14) 


[25]. 
V, V A A, B,)fA A 
g- e-e aevi as 
0 Vg VygllAy} LC DıjlAy Ay 
Using the manipulation of AC, a comprehensive Jacobian matrix with linearized equations of 
the DAE system, the state matrix As can be calculated using (15) [26]. 
As = F, — F,G,* Gy (15) 


The eigenvalues in the S-domain can be calculated since the matrix (15) has been fulfilled, which 
expresses that the system is stable if the eigenvalues of the real part are negative. 
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2.7.2. FVSI 

Fast voltage stability index (FVSI) [27] in (16) is utilized to guarantee secure bus loading since its 
value is close to 1.00. One of the buses connected to the line can cause an immediate voltage drop when the 
amount exceeds 1.00. 


20; 
FVSI,; = TAI <1 (16) 


where Z and X are the impedance and reactance of the transmission line between bus i and bus j, respectively, 
the reactive power that flows to bus j is Q;, and V; is the sending-end voltage. 


2.7.3. LQP 

The line stability factor (LQP) is a factor that guarantees a stable system when its value is less than 
1.00 [28]. At this value, the LQP ensures that no line is overloaded under all system loading conditions. The 
LQP as shown in (17). 


LQP,; = 4: (=) (SP? + 0) <1 (17) 


2.8. NSGA-II, a brief description 

In this study, the authors selected the NSGA-II algorithm [29] to solve the conflicting multi- 
objective problem as shown in (4) by finding the Pareto-optimal set. A fast non-dominant sorting approach is 
applied to this algorithm to find a group of equally good solutions closest to the Pareto optimal front. By 
calculating the crowding distance, the diversity of the Pareto-optimal solutions is maintained. Whereas 
allowing parental and offspring populations to participate in the competition increases the convergence 
nature of optimization theory, Elitism was introduced [30], [31]. Details of the procedure and critical parts of 
the NSGA-II method are given in [32]. A set of best solutions that fulfill multi-objectives, the best 
compromise solution (Best CS), is taken using a fuzzy set with full membership from the Pareto front [33]. 


3. RESULTS AND DISCUSSION 

In this study, the DFIG and UPFC are optimally installed on the grid to control the system's stability 
and security by maximizing the load bus system (Max. LBS). Investigate the Max. LBS during Min.Pioss is 
done with adequate numerous system's security and stability using the NSGA-II techniques vis: i) base case 
(without DFIG or UPFC); ii) case-1 (with DFIG only); iii) case-2 (with UPFC only); and iv) case-3 (with 
DFIG and UPFC). The effectiveness of the proposed approach for all the cases was simulated using a 
modified standard test system, IEEE 14-bus [34], [35], and the Indonesia Java-Bali 24-bus system as a 
practical test system [21], [36]. 


3.1. A modified standard test system, IEEE 14-bus 
3.1.1. Base case (without DFIG or UPFC) 

In this case, DFIG and UPFC are not connected to the grid, so all system stability is not examined. 
Then the system load is enhanced to obtain Max. LBS and Min. Pioss was 149.59% and 0.1625 p.u, 
respectively, as shown on the Pareto front in Figure 3. The minimum Pioss obtained for the best CS is slightly 
higher than the preceding result of 1.1704 p.u, but Max. LBS is much lower at 114.40%. 
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Max. LBS (f1) 
Min. Ploss (f2) 
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Figure 3. The base case’s Pareto front 
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3.1.2. Case-1 (DFIG only) 

In Case-1, the best CS is achieved by optimal placement of one DFIG on bus 8, with 49.91 MW and 
—11.56 MVAr, to carry out the best Pioss of 0.1772 p.u, but with the Max. LBS 112.24%, which is the lowest 
in Case-1, as shown in Figure 4. When placed the DFIG on bus-14, it gained the best of Max. LBS of 
157.08% and Pioss of 0.4885 p.u. All the security and stability constraints are satisfied in this case. Moreover, 
a small signal that guarantees the grid's stability on the best CS is shown in Figure 5, stated by the S field's 
negative eigenvalues. 


T r 
Non-dominated 
Best CS 

Max. LBS (f1 


Aavo 


Ploss (%) 


i i i i i i i i i 
$00 110 120 130 140 150 160 170 180 190 200 
Load Bus System (%) 


Figure 4. Case-1’s Pareto front 


Eigenvalue 


Imaginary part 


Figure 5. Case-1’s Eigenvalues 


3.1.3. Case-2 (UPFC only) 

Figure 6 detects that by installing the UPFC at the optimal location on lines 1-5, the Pareto front 
appears an increase in Max. LBS extreme reached 179.25%, with Pioss similarly a bit high at 11.02% (0.5115 
p.u.). Meanwhile, placing the UPFC on lines 9-14 with shunt and series settings of 0.412 and -13.0707 p.u., 
respectively, is achieved, Min. Pioss of 0.1911 p.u. which is the Best Poss in Case-2, even though Max. LBS 
decreased drastically only by 118.37%. However, in Case-2, all system security and stability constraints are 
met for the two conditions in Figure 7. 
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Figure 6. Case-2’s Pareto front 
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Figure 7. Case-2’s Eigenvalues 


3.2. Case-3 (DFIG and UPFC) 

For Case-3, by optimally installing DFIG on bus-4 and UPFC on lines 9-10 in the grid, it can be 
seen that there is an increase in system loading, reaching the most extreme value (Best LBS) of 193.71% with 
a high Poss of 15.31% (0.7683 p.u) as illustrated in Figure 8. The increase was achieved with DFIG capacities 
of 95.37 MW and -25.89 MVar with UPFC settings of 0.5 and -60.1915 p.u, for shunt and series, 
respectively. Intermediate results obtained on Best CS with Max. LBS and Min. Pioss is 159.75 % and 0.3261 
p.u, respectively, by optimal placement of DFIG with a size of 67.24 MW, -29.49 MVar on bus 14 and UPFC 
on lines 1-5 with shunt and series settings 0.5 and 20 p.u, respectively. Figure 9 shows that the system 
stability is fulfilled in Best CS, with all real parts of the eigenvalues being negative. While Figure 10 denoted 
the security constraints based on the FVSI and LQP for all cases in a safe condition with values less than 
1.00. 
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Figure 8. Case-3’s Pareto front 
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Figure 9. Case-3’s Eigenvalues 
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Figure 10. FVSI and LQP for all Cases 


3.3. A modified practical test Indonesia Java-Bali 24-bus system 

A simulation has been carried out on a modified Java-Bali 24-bus Indonesian practical test system as 
despicted on Figure 11 [36], [37] to prove the proposed method. Figure 12 shows the Best LBS obtained at 
165.87% with a Poss of 2,545% (1,3480 p.u.) by the optimal installation of DFIG 138.91 MW, -62.99 MVar 
on bus 13 and by placing UPFC on the exact location at lines 18-19 with shunt and series settings of 1.00 and 
-62.99 p.u., respectively. Best CS was obtained with ploss reduced by almost half, of 2.5421 p.u. But with a 
decrease in Max. LBS is not too large of 144.09% compared to the previous result, with the system stability 
is guaranteed, as shown in Figure 13 with a negative eigenvalue. The system security is also certified for all 
cases for this practical test system. It is reflected in their values less than 1.00 for the FVSI and LQP indices, 
as shown in Figure 14. 
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Figure 12. Case-3’s Pareto front of the Indonesia Java-Bali 24-bus system 
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Figure 14. FVSI and LQP of all cases for the Indonesia Java-Bali 24-bus system 


4. CONCLUSION 

A new method to improve system dynamic performance involving various security and system 
stability constraints has been successfully tested in this study with the optimal placement of DFIG and UPFC 
on the grid. The multiobjective issue was dealt with involving Max. LBS for Min. Pioss uses the NSGA-II 
technique, applied to standardized test systems and practical test systems to improve the security and stability 
of network-connected DFIG and UPFC, with all constraints satisfied. Bus voltage violation limits and line 
thermal limits are considered system safety constraints. Various system stability constraints were evaluated, 
including the voltage stability index (FVSD), the line stability factor (LQP) reaching a value of less than 1.0, 
and the small-signal stability expressed by its eigenvalues attaining negative values. The test results are 
proven by presenting the Pareto front in the best LBS (fı) condition and the best Pioss (f2) with a relatively 
extreme score of 193.71% and 0.2316%, respectively, for the standard 14-bus test system. At the same time, 
a fuzzy-based mechanism is used to obtain the best CS as an alternative to the two previously obtained 
conditions, which are 159.75% and 0.3261%, respectively. Because the type of FACTS used in this study is 
UPFC which has the most expensive investment cost among all types of FACTS devices, this research needs 
to be continued by considering additional objective functions to minimize FACTS costs. 
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